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Abstract We study possible impact of a softening of 
the equation of state by a phase transition, or appear- 
ance of hyperons, on the spin evolution of of isolated 
pulsars. Numerical simulations are performed using ex- 
act 2-D simulations in general relativity. The equation of 
state of dense matter at supranuclear densities is poorly 
known. Therefore, the accent is put on the general cor- 
relations between evolution and equation of state, and 
mathematical strictness. General conjectures referring to 
the structure of the one-parameter families of stationary 
configurations are formulated. The interplay of the back 
bending phenomenon and stability with respect to ax- 
isymmetric perturbations is described. Changes of pulsar 
parameters in a corequake following instability are dis- 
cussed, for a broad choice of phase transitions predicted 
by different theories of dense matter. The energy release 
in a corequake, at a given initial pressure, is shown to be 
independent of the angular momentum of collapsing con- 
figuration. This result holds for various types of phases 
transition, with and without metastability. We critically 
review observations of pulsars that could be relevant for 
the detection of the signatures of the phase transition in 
neutron star cores. 
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1 Introduction 

The equation of state (EOS) above a few times nuclear 
density is a main mystery of neutron stars. Our degree of 
ignorance concerning EOS of neutron star core is there- 
fore high. In particular, different theories of dense mat- 
ter predict various types of softening of the EOS above 
some density, due to phase transition to exotic phases of 
hadronic matter (meson condensates, quark matter), or 
presence of new baryons, such as hyperons (for a review, 
see Glendenning 2000). 
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Since their discovery in 1967, thanks to the excep- 
tional precision of their timing, isolated radio pulsars 
are cosmic laboratories to study physics of ultradense 
matter. Such a pulsar looses its angular momentum, J, 
due to radiation. In response to the angular momentum 
loss, J < 0, the pulsar changes its rotation frequency 
/ = 1/period = £2/2it (measured in pulsar timing) and 
increases its central density and pressure p c , P c . The re- 
sponse of the stellar structure to J depends on the equa- 
tion of state (EOS) of the neutron star core, and is sensi- 
tive to appearance of new particles or of a new phase at 
the neutron star center. In particular, crossing the phase- 
transition region by p c is reflected by specific "nonstan- 
dard" behavior of / as a function of time. In the present 
paper we analyze the relation between the pulsar spin 
down and the EOS. This problem was studied previously 
by many authors, who considered various types of EOSs 
of dense matter predicted by different theories (Weber 
& Glendenning 1991, 1992; Glendenning, Pei, & Weber 
1997, Hcisclbcrg & Hjorth- Jensen 1998, Chubarian ct al. 
2000, Spyrou & Stergioulas 2002). In the present paper 
we summarize and review results obtained using a dif- 
ferent approach. Namely, through a flexible parametriza- 
tion of dense matter EOS, we are able study the prob- 
lem of spin down evolution, including stability and the 
back bending phenomenon (spin up accompanying angu- 
lar momentum loss), for a very broad, continuous family 
of EOSs with and without phase transitions, and select 
those features which are genuine properties of the EOSs. 
On the other hand, the use of parametrized EOSs al- 
lows to perform very precise 2-D calculations, crucial for 
reliable analysis of stability, and avoiding usual severe 
limitations resulting from the lack of precision. 

The plan of the present paper is as follows. In Sect. [21 
we review effect of phase transitions on the the equation 
of state (EOS) of dense matter under the assumption of 
full thermodynamic equilibrium. Then in Sect.|3|we re- 
mind stability criteria for relativistic rotating stars. Im- 
portance of stability for the back bending phenomenon is 
discussed in Sect.QJ In Sect.[S]various cases of evolution- 
ary tracks of isolated pulsars, corresponding to different 
EOSs of dense matter, are studied, and general proper- 
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Fig. 1 Two types of the EOS with 1st order phase transition. 
Thick solid lines: stable pure N and S phases. Mixed phase 
is represented by a dash-dot line. For further explanation see 
the text. From Bejger et al. (2005) 

tics of the structure of stable families of stationary con- 
figurations are deduced. Corequakes resulting from in- 
stability of rotating configurations are discussed in Sect. 
In Sect. we show how a phase transition in neutron 
star core can affect the time evolution of the pulsar ro- 
tation period. Some remarkable features of changes in 
pulsar parameters resulting from a corequake are dis- 
cussed in Sects. |HE1 Finally Sect. 1101 contains discussion 
and conclusions. 

Results reviewed here were obtained in last few years 
within a fruitful collaboration with M. Bejger and E. 
Gourgoulhon (Zdunik et al. 2004, 2006a, 2006b, Bejger 
et al. 2005). We present these results in the most general 
terms, and summarize them in the form of conjectures 
(which in the future may become strict theorems), which 
on the one hand are general (i.e., do not depend on a 
specific dense matter model), and on the other hand seem 
to us useful for further numerical studies of the intimate 
relation between EOS and neutron star dynamics. 



2 First order phase transition in thermodynamic 
equilibrium 

Two possible cases of the EOS softening by a 1st order 
phase transition are shown schematically in Fig. ^ Sim- 
plest case is that of a first order phase transition from a 
pure N (normal) phase to a pure S (superdense) phase. It 
occurs at a well defined pressure Po, and is accompanied 
by a density jump: pn < ps- Such a type of phase transi- 
tion occurs for sufficiently strong pion or kaon condensa- 
tion, and is also characteristic for many models of quark 
deconfinement. For a sufficiently small surface tension at 
the N-S interface, a < (7 cr i t , one has to contemplate a 
more complicated pattern of the N-S phase coexistence. 
Above some P™ a mixture NS is preferred over the pure 
N phase, with fraction of the S phase increasing from 
zero to one at P = P™, and above P™ pure S phase is 



preferred. Such a situation might be possible for meson 
(kaon or pion) condensations or quark matter, provided 
g < CT cl it. Alas, our knowledge of physics of superdense 
matter is not sufficient to decide which case (pure or 
mixed) type of first order phase transition is actually 
realized in dense neutron star cores. 



3 Stability of hydrostatic equilibria 

Neutron stars are assumed to be built of an ideal fluid. 
Therefore, stresses resulting from elastic shear, electro- 
magnetic field, viscosity, etc. are neglected. 
In our case, stability criteria are: 

Stability criteria - nonrotating stars Nonrotating con- 
figuration form a one-parameter family C(x), and are 
labeled by x = p c or i = P c . One considers stabil- 
ity of these 1-D configurations with respect to radial 
perturbations and finds that a configuration is stable if 
dM/dx > and unstable if dM/dx < (Harrison et al. 
1966). 

Stability criteria - rotating stars. They were precisely 
formulated by Friedman et al. (1988). Stationary rotat- 
ing configurations form a two-parameter family: C(x, fi). 
Their stability with respect to axially symmetric pertur- 
bations can be checked using any of the following crite- 
ria, involving baryon (rest) mass Mb or stellar angular 
momentum J : 

(1) configuration is stable if {dM/dx)j =const . > and 
unstable if (dM/dx) j =CO nst. < , 

or 

(2) configuration is stable if (dJ/dx)M b =const. < and 
unstable if (dJ/dx)M b=con st. > . 



4 Back-bending and stability of rotating 
configurations 

Zdunik ct al. (2006a) have shown examples of the EOS 
with 1st order phase transitions between the pure N and 

5 phases, as well as mixed-phase transitions, for which 
one of the two situations occurs: 

— 1. all configurations are stable but the back bending 
exists. 

— 2. the phase transition results in the instability region 
for rotating stars. 

To distinguish between these two cases we can look at 
the behavior either of the curves M^(p c )j or J(/)j\/ b ; the 
quantity fixed along a sequence is indicated by the lower 
index. For such curves the instability criterion directly 
applies. However, is not so easy to detect the instability 
using the 1(f) M b plot, where / = J/J? is stellar moment 
of inertia. The unstable case may look similar to that 
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Fig. 2 Total angular momentum versus rotation frequency / 
(left panels), and moment of inertia / = J/fi versus / (right 
panels), for EOSs MSt and MUn of Zdunik et al. (2006a). The 
stability criterion is easily applied to left panels. It is clear 
that for the MSt EOS back bending feature is not associated 
with an instability, with all configurations being stable. On 
the contrary, the MUn EOS produces back bending with a 
large segment of unstable configurations. Simultaneously, the 
/(/) curves for both EOSs are very similar and apparently 
show very similar back-bending shapes. 



corresponding to a fully stable sequence. An example is 
presented in Fig. [21 where we plot two functions: J{f)M h 
and 1(f) Mh for the two models of EOS (MSt and MUn) 
with softening due a mixed-phase segment. Upper pan- 
els correspond to the stable model MSt and lower to the 
MUn model for which exist region of unstable configu- 
rations. The difference between the MSt and MUn cases 
is clearly visible in left panels (J(f)tA h curves), where 
we can easily recognize the instability region for MUn 
model, by applying the condition (dJ/dp c ) M > 0, and 
keeping in mind that p c is monotonic along this curve. 
However, the upper and lower right panels (I(f)ivi h ) are 
quite similar, without any qualitative difference. Simul- 
taneously, the back bending phenomenon in the MUn 
case is very large (almost by 200 Hz), and obviously the 
phase transition results in an instability (increase of J for 
increasing p c ). An yet, I(p c ) is a monotonic decreasing 
function all the time. 
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Fig. 3 Examples of the evolutionary tracks in the J — / plane 
for an EOS which produces back bending but does not lead 
to splitting of the tracks into disjoint stable branches, sep- 
arated by unstable segments. Curves are labeled by baryon 
mass (constant along a track) in solar masses. Thick segments 
correspond to back bending phenomenon. A small dotted ter- 
mination at the lower end of the 1.55 curve corresponds to a 
collapse into black hole after reaching the minimum J allowed 
for stationary configurations. From Zdunik et al. (2006a). 
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Fig. 4 Examples of the evolutionary tracks in the J — / 
plane for an EOS which produces unstable segments. Curves 
are labeled by baryon mass (constant along a track) in so- 
lar masses. Thick segments correspond to back bending phe- 
nomenon. From Zdunik et al. (2006a). 



5 Invariance of structure of (one-parameter) 
families {Cx} 

We select a one-parameter family (represented by a curve) 
from a two-parameter set of stationary configuration by 
fixing one of parameters, denoted by X, so that within 
this family X = const. Here X = Mb, J, /. Examples 
of such one-parameter families are shown in Fig. |21 El 



As shown by Zdunik et al. (2006a), studying stability of 
these one parameter families reveals a very interesting 
invariance property. This property can be generalized 
to include also the special case of marginally unstable 
configurations, which correspond to the inflection point 
(dM/dx) J —const . 

= and (d 2 M/dx 2 ) J=const . = 0. For 
the Mb = const, families the marginal instability corre- 
sponds to (dJ/dx) Mb =const. = and (d 2 J/dx 2 ) Mb =const. 
I). 
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Only stable configurations are astrophysically rel- 
evant. In what follows, we use the term "unstable seg- 
ment" of one parameter family in the restricted sense: 
we will always mean a segment composed of unstable 
configurations which is bounded on both sides by stable 
configurations. A remarkably general invariance property 
can be formulated as three conjectures: 

(1) All stable {C} sta t remains all stable {Cx} ro t 

(2) If {C} s tat contains unstable segment then every {Cx}n 
contains unstable segment too 

(3) If {C} s tat contains a marginally unstable C then each 
{Cx}mt contains a marginally unstable Cx 

These conjecture are based on very precise numer- 
ical 2-D simulations performed for hundreds of EOSs 
with softening due to phase transitions. They can also 
be rephrased in a more compact form: 

1' single family of stable static configurations <^> single 
family of stable rotating configurations 

2' two disjoint families of stable static configurations sep- 
arated by a family of unstable configurations O two dis- 
joint families of stable rotating configurations separated 
by a family of unstable configurations (constant Mb, or 
constant J, or constant /) 

Therefore, topology of the one-parameter families 
of hydrostatic stationary equilibrium configurations is a 
genuine feature of the EOS, not altered by rigid rotation. 
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Fig. 5 Evolutionary tracks of an isolated pulsar loosing J 
in the J — f plane, when a softening of the EOS due to a 
phase transition imply instability regions (dotted segments). 
Arrows lead from unstable configuration to a stable collapsed 
one, with the same baryon mass and angular momentum. 
From Zdunik et al. (2006a). 



6 Instabilities and starquakes 

Entering the unstable segment of the evolutionary track, 
during pulsar spin down, leads to a discontinuous transi- 
tion (starquakc) to a stationary configuration on another 
stable segment of the evolutionary track. This is illus- 
trated by examples of evolutionary tracks in the J — f 
plane in Fig. |S] It is assumed that there is no matter 
ejection during the quake and that the timescale of it 
is so short that the angular momentum loss (by radia- 
tion of electromagnetic and gravitational waves) can be 
neglected. Therefore, a spinning down pulsar reaches in- 
stability point and collapses (with spinning up!) into a 
new stable configuration C with the same baryon mass 
and angular momentum, and then continues its evolution 
moving in the J — f plane (Fig. . The unavoidability 
of a corcquake is best seen by showing the track in the 
M — J plane in Fig. The only way out from C\ is 
collapse to Cf, keeping Mb and J constant. The energy 
released during the corequakc is AE = (M ; — M f )c 2 . Let 
us stress, that the corequake is accompanied by differen- 
tial rotation, redistribution of the angular momentum, 
and breaking of the axial symmetry. It is a nonstation- 
ary phenomenon which is to be studied using methods 
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Fig. 6 Gravitational mass of the star as a function of its 
angular momentum, for a fixed baryon mass of the star. Cen- 
tral density is increasing along the curve as indicated by the 
arrows. The upper (dotted) segment corresponds to the un- 
stable configurations. From Zdunik et al. (2006a). 



beyond those applied to describe 2-D rigid rotation in 
the present paper. 
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Fig. 7 The evolution of the pulsar period P when the energy 
loss is described by the magnetic dipole breaking with a — 4. 
Solid curves - results obtained for an EOS with a phase tran- 
sition, for different baryon masses. Dotted lines correspond to 
a standard model with / = 1(0) = const. The unit of time is 
1000 yrs/K28, where K28 = ft/10 28 (in the c.g^s. units), and k 
is a parameter in the right-hand-side of Eq. Q . From Zdunik 
et al. (2006a). 



7 Back bending and pulsar timing 

A spinning down isolated pulsar looses its angular mo- 
mentum and energy. The energy balance is described by 
the standard formula 

Mc 2 = -nQ a , (1) 

where the right-hand-side describes the energy loss due 
to radiation. The standard assumption is that the In- 
dependence of the moment of inertia I can be neglected, 
so that during spin down / = I(Q = 0) = const. How- 
ever, in the phase transition epoch the fl dependence of 
I is crucial for the evolution track, timing, and stability. 
Standard approximation / = I(fl = 0) = const, leads to 
a severe overestimate of the spindown rate and an under- 
estimate of pulsar age. An example visualizing this effect 
is shown in Fig.[7J We assume that a P = 10 ms pulsar is 
observed at time = 0. As we see, evolution back in time 
using / = const, (dotted lines, corresponding to different 
Mb) can be very misleading if a phase transition region 
was crossed in the past. Let us notice that in that "phase 
transition epoch" the breaking index n = Qflju 1 was 
reaching huge values. Breaking index 3> 1 is a signature 
of a phase transition taking place at the pulsar center. 



250 r 1 1 1 1 1 1 1 1 r 




0.5 1 



J [GM|/c] 

Fig. 8 Changes of stellar parameters of a rotating isolated 
neutron stars, due to a corequake which occurs after the star 
reaches an unstable configuration. An example from Zdunik 
et al. (2006a). 



8 Spin up, shrinking of radius, and energy 
release in starquakes in pulsars 

The changes in /, equatorial radius R eq , and the en- 
ergy release can be evaluated by comparing the rele- 
vant quantities in the stationary, rigidly rotating states 
(initial and final), C; — ► Cf. Assumed conditions are 
Mb.i = Mb.f, Ji = J[. The energy release is then given 
by AE = -AMc 2 . 

Examples of numerical results obtained in 2-D sim- 
ulations are shown in Fig. |H| A remarkable feature is a 
very weak dependence of AE on Jj. Therefore AE can 
be calculated using 1-D code for non-rotating stars and 
this gives excellent prediction (within better than 20%) 
even for high Ji, when the stars are strongly flattened 
by rotation! This has important practical consequence: 
there is no need for 2-D simulations to get reliable esti- 
mate estimate of the energy release, 1-D calculation for 
non-rotating stars is sufficient. 



9 Metastability and two types of corequakes 

The case of spherical, non rotating neutron stars was 
studied long time ago (Haensel et al. 1986, Zdunik et al. 
1987). Corequakes in rotating neutron stars were recently 
studied by Zdunik et al. (2006a, b). 

First order phase transition allows for metastability 
of the N phase at P > Pq (sec Fig. EJ. Depending on 



P. Haensel and J.L. Zdunik 



t 

Fcrit 
nucl 








/ 1 

/ 




Po 


/ 

/ 




N / 







Pn Pnucl Pcrit Ps P 



Fig. 9 A schematic representation, in the p — P plane, of 
an EOS with a 1st order phase transition. Solid segments: 
stable N and S phases in thermodynamic equilibrium. Dash- 
dot segment: metastable N phase. The S phase nucleates at 
■Pnucl- At Pcrit nucleation of the S phase is instantaneous, 
because the energy barrier separating the N phase from the 
S one vanishes. 



the magnitude of the density jump at P = Po between 
stable N and S phases, one has to distinguish two types 
of 1st order phase transitions. 



One finds that in all cases AE is, to a very good 
approximation, independent of J of collapsing configu- 
ration, provided the energy release is calculated at fixed 
overcompression in the center of the metastable N-star 
core, SP = (P c — Pq)/Pq. Generally, a starquake is trig- 
gered for P c = Pnucl, when the S phase nucleates, which 
initiates a transition depicted in Fig. I1UI Independence 
of AE(SP) from J of collapsing configuration is clearly 
seen in Fig. ^2 All points lie on the same curve, which 
coincides with that obtained for nonrotating case. 

Independence of the energy release of J; is valid for 
both weak, moderate, and strong 1st order phase tran- 
sition. It has a paramount practical importance: simple 
1-D calculations for non-rotating stars are sufficient to 
get correct value of the energy release for even rapidly 
rotating stars. 




Fig. 10 Transition from a one-phase configuration C with 
a metastable core of radius rN to a two-phase configuration 
C* with a S-phase core of radius rs- The two configurations 
have the same baryon number (baryon mass) and the same 
angular momentum. From Bejger et al. (2005). 
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Fig. 11 (Color online) The energy release due to a corequake 
of rotating neutron star as a function of the overpressure 8P 
at the center of the metastable N phase core, for a model 
of the EOS with a moderate 1st order phase transition. The 
points of different color correspond to different values of the 
total angular momentum of rotating star, but they all lie on 
a same line (Zdunik et al., 2006b, to be published). Notice, 
that rotation is very rapid, with kinetic/gravitational energy 
ratio up to T/W ~ 0.1. 



Weak and moderate 1st order phase transitions charac- 10 
terized by p s /p N < § + Po/p N c 2 . A starquake is then 
triggered by nucleation of a droplet of the S phase in a 
metastable core of N phase (see Fig. . U 8 

© 

Strong 1st order phase transitions characterized by /9 s /p N > — 6 
§ + Po/p N c 2 . Configurations with P c > P with small S " 
cores are then unstable even under full equilibrium condi- < 
tions (i.e., for no metastability). They collapse into con- 4 
figurations with a large S phase core. 
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10 Conclusions 

For isolated pulsar without phase transition in the center 
the evolution of the period during spin-down is smooth 
and monotonic. For stars with baryon mass below the 
maximum allowable one for static stars, M^^, the stel- 
lar angular momentum J is monotonically decreasing 
with decreasing rotation frequency /, dJ/d/ > 0. For 
supramassive stars, with fixed baryon mass Mb > M^^^. 
the derivative d J / df changes sign very close to the Kep- 
lerian (mass shedding) limit. However, in this region the 
effect of spin-up by the angular momentum loss results 
directly from effects of general relativity, and has noth- 
ing to do with the equation of state (Cook et al. 1992, 
Zdunik et al. 2004). The relevant value of / is then so 
close to the mass shedding limit, that the star is suscep- 
tible to many other instabilities and hence it is unlikely 
to be observed as a radio pulsar. 

A softening of the EOS of dense matter due to a 
phase transition or appearance of hyperons can affect 
spin evolution tracks of isolated pulsars. In particular, 
back bending epochs and instabilities followed by core- 
quakes may appear. One parameter families (sequences) 
of rigidly rotating configurations loosing angular momen- 
tum have a stability/instability structure identical with 
that of spherical non-rotating family. Energy release as- 
sociated with a corequake in rotating pulsar depends 
only on the overpressure in the center of its metastable 
core and does not depend on the rotation rate of col- 
lapsing configuration. This independence on the initial 
rotation rate holds universally, for weak, moderate, and 
strong phase transitions, and with or without metasta- 
bility. This means that the energy release in a corequake 
in a rotating star can be calculated, with a very good 
precision, in 1-D nonrotating model with the same pres- 
sure at the core center. 

It should be mentioned that there is still no direct 
observational evidence of the back bending phenomenon. 
Although there exist many isolated pulsars with a mea- 
sured decreasing period, they are located in globular clus- 
ters and the effect of P < is usually explained by the 
acceleration in the globular cluster gravitational field (for 
a recent review see Camilo & Rasio 2005). The effect of 
period clustering (which would be another consequence 
of the phase-transition softening of EOS) is not clear 
from observational point of view. There seems to be some 
evidence for such effect (e.g. Chakrabarty 2004) for ac- 
creting neutron stars, but not for isolated neutron stars 
considered in the present paper. However, for accreting 
neutron stars the effect of back bending is expected to 
be strongly suppressed due to the simultaneous increase 
of the stellar mass during accretion (Zdunik et al. 2005). 
There might be some indirect evidence for back bending 
hided in the evolutionary history of pulsars. Some evo- 
lutionary scenarios applied to the observed pulsars re- 
sult in unrealistically short initial period if one assumes 



standard magnetic dipole braking. This paradox could 
be explained by the existence of the back-bending epoch 
during pulsar's life (Spyrou & Stergioulas 2002, Zdunik 
et al. 2006). Finally, measurements of very large values 
of braking indices (Johnston & Galloway 1999) could be 
an indication of the closeness to / = state associated 
with back bending. However, this might be explained in 
a more standard way, by the behavior connected with the 
glitch phenomena (Johnston & Galloway 1999, Alpar & 
Baykal 2006). 

This work was partially supported by the Polish 
MNil grant no. 1P03D-008-27. 
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